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Secant Tree Calculus 

Dominique Foata and Guo-Niu Han 

, ! Abstract. A true Tree Calculus is being developed to make a joint study 

^— N ■ of the two statistics "eoc" (end of minimal chain) and "pom" (parent of 

p^ ' maximum leaf) on the set of secant trees. Their joint distribution restricted 

, to the set {eoc — pom < 1} is shown to satisfy two partial difference 

(^ • equation systems, to be symmetric and to be expressed in the form of 



< 



X 



an explicit three-variable generating function. 

1. Introduction 



C^ 



> 



Q I As was done in our previous paper, whose purpose was to evaluate 

^ I the distribution of a two-variable statistic defined on the set of tangent 

r^ ■ trees [FH13], we pursue the same goal with the secant trees, by using 

an appropriate Tree Calculus. The secant trees, defined in full detail in 
the sequel, are labeled binary trees; they are called secant, because the 
number of them with 2n nodes is equal to the secant number E2n, namely, 
the coefficient of u^"^ /{2n)\ in the Taylor expansion of secw: 

^ , n>0 ^ ' 

(N. " ^2 ^4 6 ^8 ^10 

w' ■ = 1 + — 1 + — 5 + —61 + —1385 + — -50521 + ■ ■ ■ 

g ■ 2! 4! 6! 8! 10! 

^ ' (see, e.g., [Ni23, p. 177-178], [Co74, p. 258-259]). The purpose of this 

paper is to calculate the joint distribution of the two statistics 'eoc" (end 
of minimal chain) and "pom" (parent of maximum leaf) on the set of 
secant trees. To achieve this and in particular derive an explicit three- 
variable generating function two partial difference equation systems are 
to be solved. The solutions of those systems are based on Tree Calculus, 
which consists of partitioning each set of secant trees into smaller subsets, 
and developing a natural algebra on those subsets for solving the difference 
equations. 

1.1. The secant trees. The trees in question form a subclass of the 
so-called increasing, binary trees. The latter trees having n nodes labeled 
1,2, ... ,n are in one-to-one correspondence with the n\ ordinary permu- 
tations of the sequence 12 ■ ■ -n. See, for instance, Viennot [Vi88, chap. 3]. 
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Their actual definition is next stated by using the traditional vocabulary 
on trees, such as node, leaf, child, root, ... In particular, when a node is 
not a leaf, it is said to be an interior node. 

Definition. For each positive integer n an increasing tree of size n is 
defined by the following axioms: 

(1) it is a labeled tree with n nodes, labeled 1, 2, . . . , n; the node labeled 1 
being the root; 

(2) each node has no child (then called a leaf), or one child, or two 
children; 

(3) the label of each node is smaller than the label of its children, if 
any; 

(4) the tree is planar and each child of each node is, either on the left 
(it is then called the left child), or on the right (the right child); moreover, 
the tree can be embedded on the Euclidean plane as follows: the root has 
coordinates (0,0), the left child (if any) (—1,1), the right child (if any) 
(1, 1), the grandchildren (if any) (-3/2, 2), (-1/2, 2), (1/2, 2), (3/2, 2), the 
greatgrandchildren (if any) (-7/4,3), (-5/4,3), ... , (7/4,3), etc. With 
this convention all the nodes have different abscissas. The node having 
the maximum abscissa is then defined in a unique manner. Call it the 
rightmost node; it is either a leaf, or a node having a left child, but no 
right child. 

Consider the orthogonal projections of those n nodes onto a horizontal 
axis. Reading the labels of those projected n nodes from left to right 
yields a permutation a = a{l)a{2) ■ ■ ■a{n) of 12 ■■■n. This projection 
defines a bijection of the set of all increasing trees with n nodes onto the 
permutation group ©„. 

The trees which correspond to permutations xiX2 ■ ■ -Xn having the 
property that xi > X2, X2 < xs, xs > x^, . . . , in an alternating way, 
usually named alternating, are called complete increasing. In an equivalent 
manner, an increasing tree with n node is said to be complete (see [Vi88], 
[KPP94], [St99]), if axioms (l)-(4) hold with the further property 

(5) every node is, either a leaf, or a node with two children, except the 
rightmost node, which has one left child, but no right child when n is even. 
This rightmost node will then be referred to as being the one- child node. 

Each complete increasing tree with n nodes is simply called secant 
(resp. tangent) whenever n is even (resp. odd). For each n > 1 let %n 
be the set of all complete increasing trees of size n. In Fig. 1.1 the five 
complete increasing trees from T4, accordingly the five secant trees, have 
been drawn, together with the projections of their node labels on the 
horizontal axis. Notice that each of those projections, when read from left 
to right, forms an alternating permutation of the sequence 12 34. Under 
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each tree have been calculated the two statistics "eoc" and "pom" defined 

below. 
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1.2. The main two statistics. Let t E In {n > 1). If a node labeled a 
has two children labeled b and c, define niina := min{6, c}; if it has one 
child 6, let min a := b. The minimal chain of t is defined to be the sequence 
ai — 7- a2 — )■ tta — 7- ■ ■ ■ — 7- aj-i — )> a^, with the following properties: (i) 
ai = 1 is the label of the root; (ii) for i = 1, 2, . . . , j — 2 the (i + 1)- 
st term ai+i is the label of an interior node and ai+i = minai; (iii) aj 
is the node of a leaf. Define the "end of the minimal chain" of t to be 
eoc(t) := aj. If the leaf with the maximum label n is incident to a node 
labeled k, define its "parent of the maximum leaf" to be pom(t) := k. See 
Fig. 1.1. 

Those two statistics have been introduced by Poupard [Po88] on tangent 
trees (n even). She proved that "eoc" and "1+pom" were equidistributed 
on each set %2n+i of tangent trees, also that their common univariable 
distribution satisfied a finite difference equation system; she further calcu- 
lated their generating function. In [FH13] it was proved that the equidis- 
tribution actually holds on every set T^ (n > 1) by constructing an explict 
bijection (j) of T„ onto itself with the property that: l-l-pom(t) = eoc(/)(t) 
for all t. 

Working with alternating permutations it was shown in [FH12] that 
the distribution of each statistic "eoc," "pom" on secant trees satisfied 
the same finite difference system, as introduced by Poupard for tangent 
trees, but the initial conditions were different. Their generating function 
could also be calculated. 

The next step was to study the joint distribution of the pair (eoc, pom) 
on each set 2^^, that is to say, letting 



(1.2) 



/„(m, k) := #{t G T^:eoc(t) = m and pom(t) = k} 



see whether those numbers are solutions of a partial finite difference 
equation system, and try to calculate their generating function. This 
program has been achieved in our paper [FH13] for the sets of tangent 
trees T2n+i- The purpose of this paper is then to pursue this program for 
secant trees. 
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1.3. The joint distribution. The first values of tfie matrices M2n '■= 
(/2n(?T^7 k)) {2 < m < 2n; 1 < k < 2n — 1) are listed in Table 1.1. For each 
matrix have been evaluated the row sums /2n("^7») := Ylk f'^nij^ik) = 
#{t e T2n :eoc(t) = m} (resp. column sums /2n(v k) := '}2m hnim, k) = 
#{t G %2n '■ pom(t) = k}) on the rightmost column (resp. the bottom 
row). In the South-East-corner is written the total sum /2n(», •) := 
Ylimk f'^-n-i'^^k), equal to the secant number E2n- Keeping in mind that 
alternating permutations are equidistributed with complete increasing 
trees, the identity 

(1-3) f2n{',-) = ^f2n{m,k) =E2n 

m,k 

is a consequence of the old result derived by Desire Andre [Anl979, 
Anl881]. Also, as "1+pom" and "eoc" are equidistributed, we have: 

(1.4) /2n(.,fc-l)=/2n(A:,.) {2<k<2n). 

On each entry /2n("^, ^) may be defined two partial differences with 
respect to m and k as follows: 

(1.5) A /2n(w, k) := f2n{rn + 1, /c) - /2n(w, k); 

m 

(1.6) A /2n(w, k) := f2n{rn, k + 1) - f2nim, k). 

k 

By convention f2n{'m,k) := if {m,k) ^ [2,2n] x [l,2n — 1]. Our main 
results are the following. 

Theorem 1.1. The finite difference equation systems hold: 

(Rl) A'^f2n{m,k)+4f2n-2{m,k-2) = {2 < m < k - 3 < k <2n - 1); 

m 

{R2) AV2n(TO, /c) + 4 /2n-2(m, /c) = (2 < m < /c - 1< /c < 2?i - 3). 

k 

The first two top rows and rightmost two columns of the upper triangles 
{f2n{rn, k) : 2 < m < k < 2n — 1} (n > 2) can be evaluated in function of 
the row or column sums /2n-2(»7 k), f2n-2ijni •), as is now stated. 

Theorem 1.2. We have: 

(1.7) /2„(2,A;) = /2„-2(.,A;-2) = /2„_2(A;-l,.) (3 < A; < 2?i - 1); 

[First top row] 
/2n(3, A;) = 3 /2n(2, k) {4 < k < 2n — 1); [Second top row] 

/2„(m, 2n - 1) = /2n-2(TO, .) = /2n-2(», m - 1) (2 < m < 2n - 2); 

[Rightmost column] 
/2„(m, 2n - 2) = 3 f2nim, 2n - 1) (2 < m < 2?i - 3). 

[Next to rightmost column] 
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Table 1.1: the matrices M2n {^ ^ n < 5) 
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Proposition 1.3. We further have: 

(1.8) /2(., 1) = 1; /2n(., 1) = /2n-2(., •) = £^2n-2 {u > 2); 

(1.9) /2n(., 2) = 3/2.-2(., .) = ^E2n-2 (^ > 2); 

(i?3) A2/2n(m,.)+4/2n-2(m,.) = (2<m<2n-2); 
(i?4) A2/2n(.,A;)+4/2„_2(.,A:) = (1 < A; < 2n - 3). 

k 

Theorem 1.4. Tiie previous two theorems and Proposition 1.3 provide 
an exphcit algorithm for calculating the entries of the upper triangles of 
the matrices M2n = (/2n(^7 k)) (n > 1). 

The entries of the lower triangles {f2nijn,k) : 1 < k < m < 2n} 
in Tables 1.1 have been calculated directly by means of formula (1.2). 
Contrary to the upper triangles we do not have any exphcit numerical 
algorithm to get them; only the entries situated on the three sides of those 
lower triangles can be directly evaluated, as shown in Section 6. See, in 
particular. Table 6.1. For the upper triangles, we can derive an explicit 
generating function, as stated in the next theorem. 

Theorem 1.5. The triple exponential generating function for the upper 
triangles of the matrices (/2n("^, k)) is given by 

rp2n — k — l „ k—m — 1 ^m—2 

(1.10) V f2n{m,k)— — — — — — 

■^^ I2n — k — 1)1 ik — m — In im — 2)\ 

2<m<k<2n-l ^ ^ ^ ) \ J 

cos(2y) + 2 cos(2(a; — z)) — cos{2{z + x)) 
2 cos^(x + y + z) 

The right-hand side of (1.10) is symmetric in x,z. Hence, the change 
X -^ z iiv the left-hand side of (1.10) shows that 

(1.11) /2n(2n + 1 - /e, 2n + 1 - m) = /2n(m, k). 

The upper triangles already mentioned are then symmetric with respect 
to their counter-diagonals. This result can also be extended as follows. 

Theorem 1.6. Let Up(2?i) he the set of all (m, k) from [2, 2n] x [1, 2n — 1] 
such that either m — l<k, (to, k) = (3, 1), or {m, k) = (2n, 2n — 2). Then 
(1.11) holds for every (to, k) G Up(2?i). 

Those theorems and proposition are proved in the next sections, once 
the main ingredients on Tree Calculus have been developed, as done in 
the next section. 
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2. Tree Calculus 



We adopt the following notations and conventions: for each triple 
{n,m,k) let T2n,m,fc (resp. %2n,m,,, resp. %2n,,,k) denote the subset of 
%2n of all trees t such that eoc(t) = m and pom(t) = k (resp. eoc(t) = m, 
resp. pom(t) = k). Also, symbols representing families of trees will also 
designate their cardinalities. With this convention '^n,m,k '■= i^^n,m,k- 
The matrix of the upper triangles ^2n,m,k {2 < m < k < 2n — 1) will be 
denoted by Upper (2^27,,). 

Subtrees (possibly empty) are indicated by the symbols "O?" "V"? or 
"D." The notation "]^" (resp. "0," resp. "^") is used to indicate that 
the subtree "□" (resp. "O?" resp. "V") contains the one-child node or 
is empty. Letters occurring below or next to subtrees are labels of their 
roots. The end of the minimal chain in each tree is represented by a bullet 

In the sequel certain families of secant trees will be represented by 
symbols, called trunks. For example, the symbol 




A- 

is the trunk that designates the family of all trees t from the underlying 
set l2n having a node 6, parent of both a subtree "□" and the leaf m, 
which is also the end of the minimal chain. Notice that unlike the secant 
trees which are ordered, the trunk is unordered, so that 





(2.1) 

When the subtree "Q" contains the one-child node or is empty, we let 





Am ■■ 

be the family of all trees t from the underlying set T2n having a node 6, 
parent of the right child "1X[" and the left child leaf m. Let A{-k) be the 
set of all trees from A such that "Q" is not empty and does not contain 
the one-child node. We have the following decomposition 

(2.2) A = A(v.) + A(^). 
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When a further condition (C) is imposed on a trunk A we shaU use the 
notation [A, (C)]. For example, the symbol 



B 



is the trunk that has the same characteristic of the trunk A as above 
with the further property that the node labeled a belongs, neither to the 
subtree of root 6, nor to the path going from root 1 to b. 

In our Tree Calculus we shall mostly compare the cardinalities of certain 
pairs of trunks, as shown in the following two examples. 




Example 1. The two trunks 



771+1 




o 

'm+2 




,.0 

m+l 



m] and [ I , m + 2] 

5, by using the 
Example 2. To compare the cardinalities of the following two trunks 



have the same cardinalities, by using the bijection ( ^^ ^^ ^^^ ) . 




we decompose them according to the location of the one-child node: 

, 2n VV , 2n Vv , 2n vv , 2n 




:=CW + C(^)+C(^) + C(0); 




D= P + 

In the above seven trunks of the two decompositions the symbols "CH," 
"V," "O" without slash are not empty and do not contain the one-child 
node. Furthermore, note that D — D{-k) + D{^) + -D(G)) s-^so holds. 



8 
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Pivoting the two subtrees on each of the nodes k, (fc + 1), (/c + 2) in C{-k) 
(resp. in -D(*)) yields the same trunk, as in (2.1). We then say that subtree 
pivoting is permitted on those nodes. Furthermore, the three subtrees "O7" 
"I I" and "V" in (^{-k) play a symmetric role, while in D{'k) only "O" and 
"V," on the one hand, and "Q" and "V," on the other hand, have a 
symmetric role. Hence, 

CW = 2 £>(•). 

Now, remember that in each secant tree the one-child node lies at 
the rightmost position. This implies that subtree pivoting on each of the 
ancestors of the one-child is not permitted. Thus, subtree pivoting in C(^) 
and -D(^) is only permitted on nodes k + 2 and k + 1. On the other hand, 
the two subtrees "O" and '"O" play a symmetric role in C(^), but not 
in i:>(^). Hence, 

C7(^)=2D(^). 

In C(lX[) and -D(l\[) the two subtrees "O" and "v" play a symmetric role. 
Moreover, subtree pivoting on nodes k + 2 is permitted in C(]^), but not 
inD(^). Hence, 

so that 

C-2D = Cm. 

3. Proof that (Rl) holds 
The decomposition 





2^2n,m,fc = I + [ I , m + 1] 

means that in each tree from %2n,m,k the node {m + 1) is, or is not, the 
sibling of the leaf m. In the next decomposition the node m is, or is not, 
the parent of the leaf (to + 1): 




'^2n,m+l,k — I + [ I jfTi]. 

As already explained in Example 1, Section 2, we can write: 



^ ^2n,m,k ^2n,m-\-l,k ^2n,m,k ~ 
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SO that 



2n,m,k 



{^2n,m+2,k ~ '^2n,m+l,k) — {^2n,m+l,k ~ '^2n,m,k) 



m+2 




m+1 



m+1 




o 

m+2 rn+1 




m \ /m+1 



B 



C 



+ 
+ 




D 



Depending on the mutual positions of nodes to, (m + 1) and (m + 2) the 
further decompositions prevail, as again k still remains attached to (2n): 

A^ r + [ I ,m]:=Ai+A2; 

m] := Si + S2; 

m+f\>^P 

+ [ I ,m + 2] :=Ci + C2 + C3; 




m + 2] ■=Di + D2. 

In the above decompositions the subsets Si and Ci are identical. 
Furthermore, A2 = C3 and i?2 = -D2- Accordingly, A^l2n+i,mfc = 

m ' ' 

Ai—2 B1—C2+D1. A further decomposition of those four terms, depending 
upon the occurrence of the one-child node, is to be worked out: 

Ai:=AiW + Ai(H) + Ai(©); 
Si:=SiW + Si(®); 

C2:=C2W + C2(H) + C2(®); 
Si:=SiW+Si(^) + Di(0). 

Now, -Bi(*) can be decomposed into 

siW= r = r + T'" 

= Si,iW+Si,2W, 
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where "□" and "O" are supposed to be non-empty in i?i^2(*)- By the Tree 
Calculus techniques developed in Section 2, we have Ai{-k) = 2 i?i 2(*). On 
the other hand, i?i(0) can also be written 



m+l 



Bim 




m+2 



SO that Ai(^) = 2Si(0) and Ai(Q) = Si(0). Furthermore, C2(*) = 
^iW, C2{U) = 2I)i(S) and C2(©) = I5i(©). Altogether, Ai - 
2Si -C2 + D1 = (AiW + Ai(^) + Ai(0)) - 2(Si,iW + Si,2W + 

5i(®)) - (C2W + C2(H) + C2(®)) + (I)i(*) + I5i(H) + ^i(®)) = 
Ai(©)-2Si,iW-I)i(H). 

As m is supposed to be at least equal to 2, we can write 




which shows that those two families are equal. Thus 




m+l 




m+2 



A T2n,m,A: — — 2i?i^i(*) — —2 
m 

This expression is also equal to — 4'X2n-2,m,fc-25 because in each tree t 
from Bi the nodes {m + 1) and (m + 2) are both leaves. Remove them, 
as well as the two edges going out of to, and subtract 2 from all the 
remaining nodes greater than (m + 2). The tree thereby derived belongs 

to T2n-2,m,A;-2- D 

4. Tree Calculus for proving that (-R2) holds 

With n > 3 and 2 < m < k — 1 < k < 2n — 3 we have: 



fc+i 



T. 



2n,m,k 




m] + [ 




m, /c + 1 



meaning that each tree from %2n,m,k has one of the two forms: either /c + 1 
is incident to /c, or not, and the leaf m is the end of the minimal chain. 
Using the same dichotomy. 



T. 



2n,Tn,fc+l 




m] + [ 




TO, /c]. 
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As the second terms of the above two equations are in one-to-one 
correspondence by the transposition (/c, k + 1) we have: 



In the same manner, 



^2n,m,k+2 — '^2n,m,k+l — [ 




m] := B - A. 



m\:=D- C. 



Thus, A %2n,m,k — ['^2n,m,k+2 — '^2n,m,k+l) — ['^2n,m,k+l — '^2n,m,k) 
k 

= D-C-B + A. 

The further decompositions of the components of the previous sum 
depend on the mutual positions of the nodes /c, (/c + 1), (/c + 2); 

,2n 



m 




m\ 



U 2n V ^""X/ 

k+l\/k+2 fc+K^'=+2 

+ [ Y'=,m]+ [ r ,m]+ [ \k ,m] 



A 



rV; 



fc+1 





, m\ = [ I , m,k + 2] + [ 

:=Ai + A2, 
where "O" and "□" cannot be both empty in S3 and S5. 



m\ 
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Obviously, Di = Bi and Ci = Ai, so that the sum D — C — B + A may 
be written (S2 - S4) - (C2 - A2) + {D2 - S3 - S5) - 2 B2, showing that 
the components fall in three categories: (1) B2 and S4 having one subtree 
to determine; (2) C2 and A2 having two such subtrees; (3) D2, -B3 and S5 
having three of them. 

(1) With the same notation as in (2.2) write: B2 = -B2(*) + -62(1X1) 
and B4 = S4(*) + S4(^). First, S2(*) = B^{^), but S2(H) = 2S4(^). 
Hence, 



k+2 



Bo — Ba 



54(S) 




A2(*)+A2(S) + 

7^2(0). However, 



(2) Next, C2 = C2(*) + C2(H) + C'2(®) and A2 
^2(0). As before, C2(*) = A2(*) and also C2(®) 
C2(H) = 2A2(H)- Altogether, 

C2-A2 = A2{n) = [ I' ,w]- 

(3) The calculation of the third sum requires the following decomposi- 
tions: 



D2 
B3 
B, 



D2{^) + D2(^) + D2m + ^2(®); 

Bsi^) + Bs{^) + BsiU); 

S5(*) + S5(^) + S5(H); 



By the Tree Calculus techniques developed in Section 2, especially Exam- 
ple 2, we have 



D2{*) = 2 Bs{*) = 2 B^i^y, 
D2m = 2Bsm = 2Bs{^); 
D2{U) = ^Bs{U) = ^Bs{U)- 
Hence, 



D2-Bs-B, = D2m = [ 



Now 



2n\/fc+l 

A2in) = [ r ,m] = [ 




m\ 
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where the two subtrees "O" and "v" are non-empty. The last identity 
can be rewritten: ^2(1X1) = B4{^) + £'2(®)- Hence, A^l2n,7n,k = D - 

C - B + A = (S2 - S4) - (C2 - A2) + {D2 - S3 - S5) - 2S2 

fc+2^ 




fc+l 

54(H) - MH) + ^2(®) - 2 S2 = -2 S2 = -2 [ 1'= , m]. The last 

term is also equal to — 4l2n-2,m,A:, because in each tree t from B2 the 
nodes {k + 2) and (2n) are both leaves. Remove them, as well as the two 
edges going out of (/c + 1), change (fc + l) into (2n — 2) and subtract 2 from 
all the remaining nodes greater than [k + 2). The tree thereby derived 

belongs to T2n-2,m,fc- D 

5. Proofs of Theorem 1.2, Proposition 1.3 and Theorem 1.4 

Proof of Theorem 1.2. 

First top row. The trees t from %2n such that eoc(t) = 2 contain the 
edge 2 — > 1. Remove it and change each remaining node label j by j — 2. 
We get a tree t' from %2n-2 such that pom(t') = k — 2. The second equality 
is a consequence of (1.4). [] 

Second top row. As /c > 4, the subtrees "Q" and "O" below are non- 
empty, so that 



/2n(3, /c) 





1 





/2n(2,fc)= ^ :=S= ^ = S(H 



since S(*) is empty. Therefore, A(^) = 2S(H), ^(®) = 5(H) and 
then/2,„(3,fc)=3/2„(2,A;). D 

Rightmost column. The trees t from %2n such that eoc(t) = m and 
pom(t) = 2n — 1 contain the rightmost path 2n — )■ (2n — 1) — )■ . Remove 
it. What is left is a tree t' from 'X2n-2 such that eoc(t') = m. The second 
equality is a consequence of (1.4). [] 

TVext to rightmost column. If t G l2n and pom(t) = 2n — 1, then t con- 
tains the rightmost path 2n — > (2n — 1) — >, as already noted. On the other 
hand, the node with label (2n — 2) is necessarily a leaf. The transposition 
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{2n — 2,2n — 1) transforms t into a tree ti such that poni(ti) = 2n — 2. 
Also, removing the path 2n — )■ (2n — 1) and rooting either the subtree 
2n-iK^'2n^ or the Subtree 2nr^2n-i^ Q^^^ ^j^g igg^f labeled (2n - 2) gives 

rise to two trees ^2, ta such that pom(t2) = pom(t3) = 2n — 2. [] 

Proof of Proposition 1.3. For (1.8) and (1.9) we only have to reproduce 
the proofs made in Theorem 1.2 for the first and second top rows, the 
parent of the maximum node playing no role. 

To obtain (i?4) simply write (-R2) for m = 2 and rewrite it using the 
first identity of Theorem 1.2 dealing with the "first top row." Next, (-R3) 
is deduced from (-R4) by using identity (1.4). [] 

Proof of Theorem 1.4. By induction, the row and column sums 
/2n(^7 •) and /2n(»7 k) of the matrices M2„ can be calculated by means of 
relations (1.8), (1.9) and, either (-R3), or (-R4). The first and second top 
rows (resp. rightmost and next to rightmost columns) of the matrix M2n 
are known by Theorem 1.2. It then suffices to apply, either rule {R 1), from 
top to bottom, or rule (-R2) from right to left to obtain the remaining en- 
tries of the upper triangles of M2n- D 

6. The lower triangles 

Observe that for 2n > 6 the non-zero entries of the first row, first 
column and last column of M2n are identical and they differ from the 
entries in the bottom row. For instance, the sequence 5, 15, 21, 15, 5 in Mg 
appears three times, but the bottom row reads: 16, 16, 14, 10, 5. 

By Theorem 1.2 we already know that /2n(2, /c) = f2n-2{k — 1, •) = 
/2n(fc-l, 2n-l) (3 < /c < 2n-l). We also have: f2n{'m, 1) = f2n-2{rn-l, .) 
(3 < m < 2n — 1), by using this argument: each tree t from %2n satisfying 
pom(t) = 1 has its leaf node 2n incident to root 1. Just remove the 
edge 2n — > 1. Change each remaining label j to j — 1. We get a secant 
tree t' belonging to T^2n-2- The mapping t i— )> t' is bijective; moreover, 
eoc(t') = eoc(t) — 1. As a summary, 

(6.1) /2n(2,A;) = /2n(A;-l,2n-l) = /2.(A;,l) (3 < A; < 2n - 1). 

Introduce a further statistic "ent" (short hand for "Entringer") on 
each tree t from T^ (even for tangent trees) as follows: ent(t) is the 
label of the rightmost node of t. The distribution of "ent" is well-known, 
mostly associated with the model of the alternating permutations and 
traditionally called i^n^nn^'er distribution (see, e.g., [En63], [Po82], [Po87], 
[GHZIO]). We also know how to calculate the generating function of that 
distribution and build up the Entringer triangle (EntTj,(j)) (1 < j < 
n — 1; n > 2), as done in Table 6.1. 
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For instance, the leftmost 61 on the row n = 7 is the sum of all 
the entries in the previous row (including an entry on the right!); the 
second 61 is the sum of the leftmost five entries; 56, the sum of the 
leftmost four entries; 46, the sum of the leftmost three entries; 32, the 
sum of the leftmost two entries and 16 is equal to the leftmost entry. Let 
Entn(j) = #{t G 1„ : ent(t) = j}. 



j = 




1 


2 


3 


4 


5 


6 


n = 


2 
3 


1 
1 


1 












4 


2 


2 


1 










5 


5 


5 


4 


2 








6 


16 


16 


14 


10 


5 






7 


61 


61 


56 


46 


32 


16 



Table 6.1. The Entringer distribution 
Proposition 6.1. For 2 < k < 2n — 2 we have: 



(6.2) 



/2n(2n, k) = Ent2n-2{k - 1) 



Proof. Note that in a tree t G %2n such that eoc(t) = 2n and 
pom(t) = k the leaf labeled 2n is the unique son of the node labeled fc, 
which is also the rightmost node. Let (1 = ai) —)> a2 —)> as —)•■■■ —t- 
{ttj-i = /c) — 7- {ttj = 2n) be the minimal chain of t. Form a new tree t' by 
means of the following changes: 

(i) delete the path — > aj-i — > 2n; 

(ii) for i = 1, 2, . . . , J — 2 replace each node label ai of the minimal chain 
by tti+i - 1; 

(iii) replace each other node label 6 by 6 — 1. 
The label of the rightmost node oft' is then equal to the pom(t) — 1 = k — 1. 
The mapping t i— )■ t' is obvious bijective. [] 

Finally, the upper diagonal {f2n{k + 1, /c) : (1 < fc < 2n — 1)} of the 
lower triangle in each matrix M2n can also be fully evaluated. First, in an 
obvious manner, 

(6.3) /2„(2,l) = /2n(2?i,2n-l)=0. 
Second, the identities 

(6.4) /2^(3, 2) = 2 /2^(3, 1) = /2n(2n - 1, 2n - 2) = 2 /2„(2n, 2n - 2) 

(also equal to /2n-4(», •) = -£'2^-4) can be proved as follows. To each tree 
t G ^2n such that eoc(t) = 3, pom(t) = 1 there correspond two trees, 
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whose "eoc" and "pom" are equal to 2 and 3, respectively, as illustrated 
by the diagram: 





/2n(3,l)=^"^^ ^ ^^ =/2n(3,2). 

1 1 

Likewise, each tree t G %2n such that eoc(t) = 2n and pom(t) = 2n — 2 
necessarily has its rightmost four nodes equal to (2n — 1), 6, 2n, (2n — 2) in 
that order, with b being a node less than {2n — 2). In particular, the latter 
node is its rightmost (one-child) node. To such a tree there correspond two 
trees, whose "eoc" and "pom" are equal to (2n — 1), (2n — 2), respectively, 
as illustrated by the next diagram. In particular, the node b becomes the 
rightmost node of those two such trees. 



2n-l 




2n-2 




Proposition 6.2. We have the crossing equalities: 

(6.5) f2n{k -l,k) + f2n{k + 1, k)=f2n{K k - I) + f2n{k, k + 1) , 

for3<k<2n-2. 

The involved entries are located on the four bullets drawn in the 
following diagramme. 



k-1 k k + 1 



k-1 

k 
k + 1 



Proof. Let i, j he two different integers from the set {(/c — l),/c, (/c+1)}. 
Say that i and j are connected in a tree t, if the tree contains the edge i — j, 
or if i and j are brothers and of them is the end of the minimal chain of t. 
Each of the four ingredients of the previous identity is now decomposed 
into five terms, depending on whether the nodes (/c — 1), k, {k + 1) are 
connected or not, namely: no connectedness; only k, {k + 1) connected; 
(/c — 1), /c connected; (/c — 1), (/c + 1) connected; all connected. Thus, 



fc-i 



f2n{k- l,k) 




2n 




k-1 



k + l] + [ 




2n\ /k+1 

V 
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2n 




2n 



.• fc-l 



A:-l\ /k+1 




A1+A2 + AS + A4 + A5 




A:+l 



/2n(A: + l,A;) 




2n 



fc+l 



/c-l. 



^'^V^p 




fc+1 



+ [ 




2n 



■r V.;^^, 



:= Si + S2 + S3 + S4 + S5; 



/2.(A;,A;-1) = [ 





2n 
k + 1, >-' ] +[ 



+ 
k-1] 
+ 



f'^ /fc+1 2n 




,• fc-l 



2n 







fc 




2n 



2n^^/fc+l fc 

+ [ P"\A; + 1] + [ ^'^-^ 
Ci + C2 + C3 + C4; 




f2n{k,k + l) = [ 




r+' A; - 1] + [ 




Now, the following identities hold: Ai = Ci, A2 = C4, A^ = D2, A4 = 
C2, Si = Si, S3 = S4, S4 = S3, so that J2^{A, + Bi) - J2,{C, + A) = 

(S5-S5)-(C3-A5-S2). 

As before, we may write S5 = B^i^) + S5(0), S5 = S5(*) + S5(0). 
As SsW = S5(*) and S5(®) = 2S5(®), we get: 



(6.6) S5-S5=S5(©) = 

as k is supposed to be greater than 3. 
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Next, 

C3 — i?2 = [ 

fc-l\ /fc+1 fc-1^ / fc~ 

/c + 1] - [ 

2n 

fc-l\ /fc+1 fc^ 

C3 — i?2 — ^5 = 

= (£;W-FW) + (£;(©) -F(®)) 




k-i] 






c+l fc-1 

fc-l\ /fc 
(6.7) = F(©) = 

By comparing the evaluations (6.6) and (6.7) we get: (S5 — D^) — (C3 — 
A5 — B2) = 0. This completes the proof of (6.5). [] 

The entries /2n(3,2) and /2n(2n — l,2n — 2), belonging to the upper 
diagonal of the lower triangle matrix M2n, being evaluated by identity 
(6.4), and the entries of the upper triangle being known by Theorem 1.4, 
we can apply the crossing equalities (6.5), starting with /c = 3, and obtain 
all the values of the entries of that upper diagonal. Altogether, besides 
the entries of the upper triangle, all the entries lying on the border of the 
lower triangle can be calculated, as illustrated in boldface in the following 
matrix Mg. 



Ms 



k = 


1 


2 


3 


4 


5 


6 


7 


/8(m,.) 


m = 2 






5 


15 


21 


15 


5 


61 


3 


5 


10 




45 


63 


45 


15 


183 


4 


15 


35 


50 




101 


63 


21 


285 


5 


21 


54 


86 


106 




45 


15 


327 


6 


15 


46 


82 


87 


50 




5 


285 


7 


5 


22 


46 


60 


40 


10 




183 


8 




16 


16 


14 


10 


5 




61 


fsi-.k) 


61 


183 


285 


327 


285 


183 


61 


Eg = 1385 



Table 6.1: the matrix Mg, bold-faced entries analytically evaluated. 
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7. Generating functions for the /2n(^5^) 

The calculation of the generating functions for the /2n(^7 k) is similar 
to the calculation made for the tangent tree case in a our previous paper 
[FH13]. Recall the definition and some basic properties of the Poupard 
matrix. Let G = (gij) {i > 0, j > 0) be an infinite matrix with 
nonnegative integral entries. Say that G is a Poupard matrix, if for every 
^ > 0, J > the following identity holds: 

(7-1) 9i,j+2 - 2 9i+i,j+i + gi+2,j + 4 Qij = 0. 

Remark. Last coefficient is 4, not 2 as in [FH13]. 

Let G{x,y) := E 9^,J (xytl) {y^ /jl); R^{y) := E 9^,J iv' /j'-) 

i>0,j>0 j>0 

{i > 0); Cj{x) := E 9i,j (a^V'^0 (j — 0) be the exponential generating 

functions for the matrix itself, its rows and columns, respectively. Propo- 
sitions 9.1 and 9.2 in [FH13] can be rewritten as follows. 

Proposition 7.1. The following four properties are equivalent. 
(i) G = (gij) {i > 0, J > 0) is a Poupard matrix; 
(ii) R'liy) - 2 R[^^{y) + R,+2{y) + 4 i?,(y) = for all i > 0; 
(Hi) C'^{x) - 2 q+i(a:) + C,+2(^) + 4C,(x) for all j > 0; 
(iv) G{x, y) satisfies the partial differential equation: 

d^G{x,y) ^ d^G{x,y) ^ d'^G{x,y) ^ ^(jr^ \ _ q 
dx"^ dxdy dy'^ 

Proposition 7.2. Let G{x, y) be the exponential generating function for 
a Poupard matrix G. Then, 

(7.3) G{x,y) = A{x + y) cos{2 y) + B {x + y) sm{2y), 

where A[y) and B{y) are two arbitrary series. 

The entries of the upper triangles of the matrices (M2n) (see Table 1.1) 
are now recorded as entries of infinite matrices (0*^^-*) {p > 1) by 



(7.4) oj. 



(p) ._ / 0' ifi+j=p mod 2; 

f2n{m,k), ifz+j^pmod2; 



1,3 

with m := p + 1, k := p + j + 2, 2n := p + i + j + 3. Conversely, 
i := 2n — k — 1, j := k — m — 1, p := m — 1. In particular, the first 
one Q^^' = {(^i j) {i-,j > 0) contains the first rows of the upper triangles, 
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displayed as counter-diagonals. Furthermore, a counter-diagonal with zero 
entries is placed between two successive rows. 

0W = 










1 


2 


3 


' 


4 


5 


6 


7 





/^/4(2,3) 





/6(2,5) 





/8(2,7) 





/io(2,9) 


o\ 


1 





/6(2,4) 





/8(2,6) 





/io(2,8) 







2 


/6(2,3) 





/8(2,5) 





/io(2,7) 









3 





/8(2,4) 





/io(2,6) 











4 


/8(2,3) 





/io(2,5) 













5 





/io(2,4) 















6 


/io(2,3) 

















7 


\ 














/ 





1 2 


3 4 


5 6 


7 











I 


1 


5 


61 


o\ 










1 





3 


15 


183 












2 


1 


21 


285 


••• 












3 





15 


327 














4 


5 


285 

















5 





183 


















6 


61 




















7 V 



/ 



Proposition 7.3. Every matrix Vt^'P' (p > 1) is a Poupard matrix. 
Proof. Using Definition (7.4) we have 



u 



(p) 



,(p) 



(p) 



,(p) 



i,j+2 - '^^i+lj + l + ^i+2,j +^^i,^ 

== /2n+2(TO, k + 2)-2 f2n+2{m, k + 1) 

+ /2n+2(m, k)+4 f2n{rn, k) 
= A /2n+2(TO, k)+A /2n(m, k) = 0, 



by rule (i? 2). Q 

The row labeled i of Q^^' will be denoted by O^^J and the exponential 
generating function for that row by O^^j(y) = X]i>o'^j j 2^"'/^'- Also, 
fl^P\x,y) := Xli>o ^i . v^)-^'/^' "^ill be the double exponential generating 
function for the matrix O'-^-'. 

Proposition 7.4. For ail p > 1 we have: 



dy 



dy 



21 



DOMINIQUE FOATA AND GUO-NIU HAN 
Proof. For the first identity it suffices to prove Uqj = u^j , tfiat is 

/2n(m, 2n - 1) = /2n(2, 2n + 1 - TTl) . 

This is true by the symmetry property of the Poupard triangle, as proved 
in Corollary 1.3 in [FH12]. For the second identity it suffices to prove 

oj[j = 3oJqj^^, that is 

f2n{m, 2n - 2) = 3/2n(m, 2n - 1). 

This is true by Theorem 1.2. [] 

As X and y play a symmetric role in (7.3), the solution in (7.3) may 
also be written 

G{x, y) = A{x + y) cos(2 x) + B{x + y) sin(2 x), 
so that the generating function of each matrix O*^^-* is of the form 

n(p)(x,y) = A{x + y) cos(2x) + B{x + y) sin(2a:). 

Hence, n^P\x,y) = o[,^](y) = A{y). Also, 

{x=0} 

{d/dx)n^P\x,y) = {{d/dx)A{x + y)) cos{2x) + A{x + y){-2)sm{2x) 
+ {{d/dx)B{x + y)) sin(2 x) + B{x + y)2 cos(2 x) 

and 



d 



d 



—n(p\x,y) =—A{x + y) +2B{x + y) 

ox {x = 0} ox {x = 0} {x ^ 0} 

= ±A{y) + 2Biy) 

= <^(l/)- 
By Proposition 7.4 we have A{y) = fl _i ^{y) and B{y) = (O^^j(y) — 

(7.5) n^P\x, y) = cos(2 x) 0j,'2i,.(x + y) + sin(2 x) -^^^p-iJ^ + v)- 

First, make the evaluation of Vt'^^> (x, y) . The row labeled of the matrix 
fl^^^ reads: 1, 0, 1, 0, 5, 0, 61, . . . , which is the sequence of the coefficients 
of the Taylor expansion of secy. Thus, Oq iiu) — secy. Taking p = 1 in 

(7.5) we get 

Q^^' (x, y) = sec(a: + y) cos(2 x) + sec(x + y) tan(a: + y) sin(2 x) 

(7a\ cos(x - y) 

(7-6) = j7 — ■ — r. 

cos^(x + y) 
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For further use let us also calculate the partial derivative of Q^^\x, y) 
with respect of y: 

d 1 / 

7— O*^^-* ix, y) = T- ( sin(a; — y) cosfx + y) + 2 cosfx — y) sinfx + y) 

oy cos"^(x + y) V 

(7.7) = 3I (sin(2y) + 3sin(2x)). 

2 cos-^^x + y) 

Now, define 

p-i 



(7.8) n{x, y, z) := J^ Q^^\x, y)j^^—^. 



^l'^^{x,y)j^ 
p>i 

and make use of (7.5) — (7.8): 



p—i 
n{x,y,z) =cos{2x)^n'p\^,{x + y) _ 

p>i ^^ '' 

+^'"<2-)|(i:«^'i..(-+!')(^ 

d 

= cos(2a;) O^^^ (2, x + y) + sin(2a;) 7— O^^^ (2, x + y) 

cos(x + 2/-z) 
= cos(2a:) 



cos^(x + y + z) 



2 cos"^(x + y + 2) 
cos(2a;) (cos(2(a: + y)) + cos(22)) 



2 cos"^(a: + y + z 

+ sin(2a;) (sin(2(x + y)) + 3sin(2z)) 

(7.9) = ^7 fcos(2y) + 2 cos(2(x - z)) - cosi2iz + x)) 

^ ^ 2 cos3(x + y + 2) V ' ^ ^ " " 

ip) 
By definition of the w) s we get 



™2n— A; — 1 „ k—m — l ^m — 2 

(7.10) = J2 f2n{m,ky y 



{2n-k-l)\ {k-m-l)\ (m-2)! 

the latter sum over the set {3<m + l<A;<2n — 1}. We have proved 
Theorem 1.5. 
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